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ABSTRACT 

Circumstellar disks provide the material reservoir for the growth of young stars and for planet 
formation. We combine a high-level radiative transfer program with a thermal-chemical model of a 
typical T Tauri star disk to investigate the diagnostic potential of the far-infrared lines of water for 
probing disk structure. We discuss the observability of pure rotational II2O lines with the Herschel 
Space Observatory^ specifically the residual gas where water is mainly frozen out. We find that 
measuring both the line profile of the ground lio — loi ortho-Il20 transition and the ratio of this line 
to the 3i2 — 3o3 and 221 — 2i2 line can provide information on the gas phase water between 5-100 AU, 
but not on the snow line which is expected to occur at smaller radii. 

Subject headings: accretion, accretion disks - infrared: stars - planetary systems: protoplanetary disks 
- stars: formation: pre-main sequence - X-rays: stars - radiative transfer 



1. INTRODUCTION 

The radius where volatiles sublimate in a protoplan- 
etary disk plays an important role in planet formation. 
One volatile of particular interest is water, and the place 
where it su blimates or fre ezes out onto grains is called the 
snow line (|Havashilll981l ). The sublimation temperature 
for water ranges from T ^ 110 — 170 K, depen ding on the 
water vapor pressure (cf. Eraser et al. 120011 : Podolak & 
' Zucker|200J) . At radii larger than the snow line, it is pos- 
sible t o hydrate planetesimals and protoplanets. iHavashil 
(|1981h found that the snow line is now located at about 
2.7 AU from the Sun in our own Solar System. However, 
the position can be different for a T Tauri star. Its loca- 
tion depends on the evolutionary ph ase of the disk, and 
. it can vary from 0.7 to 10 AU ( Sass elov fc Lecad [2OOOI : 
■ iLecar et al.ll2006l : iGaraud fc UiimM . 

Gas phase water is expected to be abundant between 
the water-ice sublimation temperature T ~ 110 — 170 K 
and the dissociation temperature T ~ 2500 K. Its 
rich spectrum covers almost all astronomically available 
wavele ngth bands. Using the near- infrared lin es near 
2.3^lm, iCarr et all (I200I and iThi fc BikI ([2005h found 
warm H2O in protoplanetary disk s at temperatures T > 
120 K and radi i R < .4 AU, and lCarr fc Naiitel (l2008h 
and ISalvk et"aLl (|2008f ) have detected mid-infrared wa- 
ter lines with Spitzer. The far-infrared spectrometers 
aboard the Herschel space telescope will offer the oppor- 
tunity to measure the pure rotational transitions of this 
important molecule, which are characteristic of gas tem- 
peratures of T > 60 K. Because gas phase H2O requires 
warm conditions (T > 200 K) it is expected to be abun- 
dant at small radii and thus perhaps not observable with 
the Herschel spectrometers. However, temperatures ap- 
propriate to both the occurrence and detection of H2O 
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are present over a wide range of radii because stellar radi- 
ation heats the disk atmosphere. Moreover, as we show 
in this letter, several lines of H2O can be detected by 
Herschel at very low abundances close to the mid-plane 
where the bulk of the water is frozen. The detection of 
these lines will complement the previous measurements 
of the near-infrared and mid-infrared water lines. 

2. MODELS 

The calculation is done as follows: (1) We use an X-ray 
irradiated disk code to calculate the temperature struc- 
ture and molecular abundances for a generic T Tauri 
disk. (2) The results are input into a multi-zone radiative 
transfer code that includes the calculation of the excita- 
tion of the H2O Unes. (3) A ray-tracing code is used to 
obtain line fluxes and line shapes. These three steps are 
discussed in more detail below. 

2.1. The thermal and chemical structure of the disk 

The thermal-chemical structur e of the disk is calcu- 
lated with the code described bv iGlassgold et all (I2004D 
with minor corrections and updates ( Meiierink et al.l 
12008( 1. The disk is illuminated by stellar X -rays with 
a thermal spectrum with temperature Tx — 1 keV and 
luminosity Lx — 2 x 10'^*' erg s~^. In regions where the 
X-rays are strongly attenuated, the disk is ionized by 
the decay products of ^^Al at a rate C26 = 4 x 10~^^ s~^. 
The density structure is given by the generic T Tauri disk 
model (Paola d'Alessio, private communication) with ac- 
cretion rate M = lO^^M0yr~^ and stellar parameters 
M, = 0.5 M0, = 2R0, and = 4000 K. The 
disk is flared, and the density varies continuously from 
R ^ 0.028 to >500 AU. The model density does not 
include modifications such as holes, gaps, and rims sug- 
gested, for example, by the spectral energy di stribu tions 
measured with Spitzer (e.g., DuUemond et al. l2007f ).The 
adopted model has a power law grain size distribution 
with index p = —3.5 and minimum and maximum grain 
sizes of amin — 0.005 and a„iax = 1000 //m, respectively, 
resulting in a geometric mean grain size, Qg = 2.24 /im. 
The grain size ag determines the thermal coupling be- 
tween dust and gas. 
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Fig. 1. — Temperature (top), and fractional water abundance 
(bottom) vs. perpendicular column density. The dotted blue line 
(top) shows the freeze-out temperature Tf = 110 K. The colored 
tick marks at the top of the bottom panel indicate the smallest per- 
pendicular column density for each radius where water freezes out. 
No freeze-out occurs within R = I AU. Note the non-monotone 
variation of the water abundance due to gradients in the ionization 
rate, density, and temperature. 



The lD'Alessio et al.l (|1999ll2001l) dust temperature in- 
volves a balance between heating by viscous dissipation, 
FUV/optical/IR radiation, cosmic rays and cooling by 
IR radiation. Our gas temperature is determined by 
heating by penetrating X-rays, gas-grain collisions and 
cooling by atomic and molecular lines. The gas thermal 
model relies solely on X-r ay irradiation in contrast t o 
UV ir radiatio n mode l s bv iJonkh eid ct al." (2004', '2007^, 
iKamp & DuUc mondl (|2004l ). and Nomura ct al. (2007) 



(who also include X-rays), who include UV photons. 
The UV irradiation, however, play no role in the regions 
where the submillimeter H2O emission is produced. The 
numerical approach is "1-l-lD", i.e., the vertical ther- 
mal and chemical structure is calculated independently 
at 45 radial points whose spacing increases with radius. 
However, the X-ray attenuation is calculated along the 
linc-of-sight to the star. The photons can only escape in 
the vertical direction, a reasonable approximation, since 
the opacities in the radial direction are so much larger. 

The chemistry is steady state and entirely gas phase, 
and it includes only carbon and oxygen bearing species. 
It does not treat freeze-out of H2O in detail, but we ex- 
pect that freeze-out cannot be 100% efficient. Inc omplete 
freeze -out has be e n obs erved by ISO (Boonm an et al.l 
I2003f ). and iDavisI (|2007t ) has modeled freeze-out m ir- 
radiated accretion disks. In order to treat the freeze- 
out of water, we assume a fixed gas-phase water abun- 
dance, x/(H20) = 10-8 or x/(H20) = 10"!°, in re- 
gions where most of the water is expected to freeze-out 
on grains, i.e., where the gas temperature satisfies the 
condition T < Tf with Tf — 110 K, the adopted freeze- 
out temperature. Where the model gas phase abundance 



is lower, we adopt the calculated abundance. A more re- 
alistic treatment of water freeze-out would involve treat- 
ing many poorly understood physical processes beyond 
the scope of this short report. These might include the 
growth and settling of the dust grains, molecular synthe- 
sis and transformation on grain surfaces, adsorption and 
desorption, and mixing. We are developing an improved 
thermal-chemical disk model that will address many of 
these issues. The advantage of the present model is that 
it allows us to isolate the chemical effects of freeze-out 
and the effects of grain growth and settling in a simple 
albeit preliminary way. 

2.2. Water radiative transfer 

The critical densities of the lines of interest are typ- 
ically large, in the range Ucrit ~ 10® — 10^*^ cm"^. 
The fractional abundance of gaseous water becomes sig- 
nificant at perpendicular column densities Nn > 2 — 
5 X 10^^ cm~^ (Figdl), and the excitation can be sub- 
thermal, especially at large radii {R > 10 — 20 AU). Thus 
assuming local thermodynamic equilibrium is not always 
a good approximation, and a full radiation transfer cal- 
culation is required. 

We calculate the level populations of II2O using the 
multi-zone escape p robability method /33D described by 
iPoelman fc Spaaiia C2005, 2006 ). We assume Keplerian 
rotation for a stellar mass = 0.5 Mq and turbulent 
velocities Sv=0.5, and 2.0 km (where Sv is the 1/e 
half width of line profile) . In the multi-zone formalism, 
the medium is divided into a large number of grid cells, 
each with calculated values of the volumetric hydrogen 
density, the temperatures of gas and dust, and the II2O 
abundance. The connection of every grid cell with all the 
others enables us to solve locally the equations of statis- 
tical equilibrium taking into account global information 
for the entire disk. Radiative and coUisional excitation 
and de-excitation by electrons and atomic and molecu- 
lar h ydrogen, and dust pumping are included, using data 
from lFaure fc JosselinI ([2008) and the Leiden Atomic and 
Molecular Database (LAMDA, Scho ier et al I2005D We 
adopt the dust opacities from D'Ales sio et al.l ( 20011 ). Al- 
though the radiative transfer code can deal with three 
dimensions, we only consider ID slabs as in the thermal- 
chemical calculation. 

2.3. Ray tracing 

We compute line profiles for inclinations ranging 
from 0° to 90° with the sky brightness distribution 
program SKY, which is part of the RATRAN code 
(iHogcrhcii de fc van der TakI |2000|) . We use the 2- 
dimensional cylindrically-symmetric version of SKY, and 
the resulting level populations, densities and molecular 
abundances are combined and rebinned into a single grid 
at each radius. The dynamical ranges of the physical 
and chemical properties are large in the vertical direc- 
tion, and we use 400 cells in this direction and decrease 
the cell size toward the disk midplane. Fewer grid cells 
are needed in the radial direction since the changes are 
more gradual, although smaller spacings are also used at 
small radii. 

3. RESULTS & DISCUSSION 

Fig. [1] shows the gas temperature and H2O abun- 
dance for a range of radii from i? = 1 — 40 AU plot- 
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ted versus increasing perpendicu lar column dens it y A^h . 
As explained in mo r e det ail in iGlassgold et alJ (|2004l ) 
and iMeiierink et al.l (|2008D . the stellar X-rays heat the 
gas at the top of the atmosphere as high as T ~ 
3000 - 4000 K for i? < 20 AU. Attenuation of the X- 
rays for A'h > 10^^ cm~^ produces a drop in temper- 
ature and induces transitions from atoms (or ions) to 
molecules. For example, the H to H2 transition occurs 
for A'h = 10^^ — lO^'^ cm~^. The water abundance be- 
fore the H to H2 transition is small due to the efhcient 
destruction by X-ray ionized species such as H+ and Hc+. 
Congruent with the H to H2 transition, the water abun- 
dance increases steeply. 

According to Fig. [2 freeze-out occurs for R > 1 AU 
and beyond A'h ~ 10^^ — 10^^ cm~^ depending on the 
radius (as indicated by the colored tick marks for each 
radius). The volumetric densities here are quite high, 
riH = 10^ — 10^^ cm"-^ and, despite the reduced dust 
surface area in the thermal-chemical model, the gas is 
thermally coupled to the dust. At even higher column 
densities than are shown in Fig. [T] (closer to the mid- 
plane), the temperature increases slightly due to viscous 
accretion heating. Temperatures greater than T ^ 200 K 
are needed for efficient gas phase synthesis of water by 
neutral radical reactions, and they are only attained at 
small radii R < 0.5 AU, where xh,o > lO^"' for Nh > 

1022 _ 1023 cni-2. 

We have calculated line profiles and fluxes for the low- 
excitation H2O transitions that are in the spectral range 
of Herschel Photodetector Array Camera and Spectrom- 
eter (PACS) and Heterodyne Instrument on the Far- 
Infrared (HIFI) . The results indicate that a number lines 
in the frequency bands of HIFI are detectable. The fre- 
quency bands of PACS are less favorable, and the follow- 
ing discussion is based on HIFI. The integrated fluxes 
are summarized in Tabic [T] for sources at an assumed 
distance of 140 pc. During the first Herschel observa- 
tion cycle, only the ortho lio — loi and para In — Oqo 
ground level transitions will be observed in the Herschel 
key program "Water in Star-forming Regions With Her- 
scheF (http://www.strw.leidenuniv.nl/WISH/), and we 
discuss these lines first. 



TABLE 1 

Integrated intensities / = jTMsrff [mK km s"-'-]" 



3.1. The lio — loi and In — Oqo H2O lines 
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Fig. 2. — The lio — loi ortho-H20 lines with the d'Alessio 
opacities for inclinations of 0° (left), 30° (middle), and 60° (right); 
no freeze-out (black), with freeze-out below 110 K and abundance 
a;/(H20) = 10"** (green) and x^(H20) = lO"!" (red) in the freeze- 
out zone; turbulent velocity 5v = 0.5(top), and 2.0 km s"-"^ (bot- 
tom). 

The line profiles of the 556.936 GHz lio — Iqi ortho 
and the 1113.34 GHz In — Oqo para transitions are con- 
volved with the 38" and IQ" Herschel beam sizes at these 
frequencies. Fig. [2] shows how the profile of the lio — loi 
ortho transition is affected by changes in inclination an- 
gle i, turbulent velocity 5v, and the effects of incom- 
plete freeze-out. Variations in i are shown horizontally 
across the figure, and changes in 5v vertically; the three 
lines represent no freeze-out and freeze-out that leaves 
residual water abundances of 10~^ and 10~^°. The in- 
tegrated intensities for the models shown vary between 
14 and 66 niK km s^^, and are mainly sensitive to the 
line width. The lio — loi ortho- and In — Ooo para-H20 
lines have similar integrated intensities, with line ratios 
lio ~ loi/lii ~ Ooo ~ 0.6. The small variations in this 
ratio, caused by high opacities in the lines, makes these 
lines unfavorable as diagnostics. 

Due to the rotation of the disk, the lines are broadened 
with increasing inclination, but at increasingly smaller 
integrated intensities because of the reduction of the ap- 
parent surface area of the disk. For a completely edge on 
disk, absorption by the cold outer regions {R > 100 AU) 
is important, but this effect has not been included in the 
present calculations. At an inclination of 60°, the ro- 
tational broadening extends to 5 km (best seen at 
Sv — 0.5 km s^^), which corresponds to a rotational ve- 
locity of Vr = 10 km s~^ or a radius of i? = 10 AU. The 
dominant contributions to the lig — Igi line flux originate 
therefore from i? > 10 to 100 AU. Thus this line provides 
no information on the snow line, which is located close 
to the star, near R ^ 1 AU in this model (see Fig. [T]). 

At an inclination i = 0° (face-on disk), the line width 
is determined by the turbulent velocity 6v and the line 
opacity. The center of the line is strongly self-absorbed 
for no freeze-out, since the emission produced at A'h 
1022 _ io23 cni~2 is absorbed by the upper layers. These 
regions are sub-thermally excited because the densities 
are low (nn < 10* cm^"^). Self-absorption occurs, but 
to a lesser extent, when freeze-out is included and the 
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residual water vapor abundance is (H2O) = 10^^. Even 
for a;/(H2O)=10^^", the line is still optically thick. For 
finite inclination angles, i = 30° and 60°, flux loss at 
line center can also occur for the freeze-out model with a 
water vapor abundance a;/(H20) = 10^^°. This is most 
clearly seen for the case Sv = 0.5 km s~^. This is not 
a signature of the snow line, since the flux loss occurs 
for i — 60°, for example, for velocities v < 2 km s~^, 
corresponding to radii i? > 50 AU. 



i=o° ;=30° :=60° 
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Fig. 3. — The lio — loi ortho-H20 line profile with a turbulent 
velocity of 2.0 km for inclinations (from left to right) i = 0°, 
30°, 60°, with freeze-out and a;j(H2 0) = 10"*. The model results 
are in black. For the curves in green, a Itr noise level of 0.5 mK has 
been added. The velocity resolution is assumed to be 1 km s~^. 

In Fig. [21 we test how well the shape of the lio — Iqi 
line can be measured by Herschel. The model is shown 
in black, and a la noise level of 0.5 mK has been 
added in green. The sensitivity of the Herschel detec- 
tors at 1113.34 GHz is about four times less than at 
556.936 GHz. According to the time estimator for the 
Herschel Space Observatory (http://www.esac.esa.int), 
a la noise level of 0.5 and 2.0 mK can be obtained at the 
frequencies of the lio — loi ortho and In — Oqo para tran- 
sition, respectively, in about 4 hours of observing time, 
when frequency switching is used and a spectral resolu- 
tion of 1 km s~^ is requested. Thus one should be able 
to determine the width, and to some degree the shape 
of the lio — loi ortho-H20 line. However, the effect of 
self-absorption due to high opacities will be hard to dis- 
tinguish for low turbulent velocities {5v < 1.0 km s~^) 
at this resolution. The para line is not a good candidate 
as the ortho transition for obtaining a line shape at the 
same noise level because of the decrease in sensitivity. 

3.2. Line ratios 

The line profiles and integrated intensities of the lio — 
loi and its ratio with the In — Ooo transition are af- 
fected by inclination angle, turbulent broadening, and 
the residual water vapor abundance after freeze-out. The 
dust opacity, which is influenced by grain growth and set- 
tling, can also play a role, although our study of these 
effects indicates that this is less important than the vari- 
ables just mentioned. Thus the models are degenerate 
in the sense that different parameter sets give the same 
line shapes and intensities, due to the high opacities in 
these lines. Unfortunately, these are the only lines to be 
observed in the WISH program. Observing higher exci- 
tation transitions would help, but the decrease in sen- 
sitivity with frequency means that only line fluxes and 
not line shapes can be determined by Herschel observa- 
tions. At a velocity resolution of 1 km s~^, it would take 
Herschel about 4 hours to get the la noise level down to 
2 mK for transition frequencies around ~ 1100 GHz. In 
a similar observation around ~ 1700 GHz, the noise level 



would be 8 mK. 
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Fig. 4.— The 3i2 - 3o3/liO - loi versus the 221 - 212/I10 - loi 
ratio for models without freeze-out (black), models with gas-phase 
abundances in the freeze-out zone of xy(H20) = 10~* (green) and 
IQ—iO (orange). The different lines are for inclinations i = 0° 
(solid), 30° (dotted), and 60° (dashed). 

Inspection of Table 1 shows that higher transitions may 
be observable with Herschel. We will now focus on the 
3i2 — 3o3 transition at 1097.37 GHz and the 221 — 2i2 
transition at 1661.01 GHz, although other lines may also 
have diagnostic value. Fig. [4]plots the ratios of these two 
lines to the fundamental lio — loi transition, and shows 
that the ratios can distinguish between different levels of 
water vapor in the freeze-out zone. The model results 
for no freeze-out are found in the right upper part of the 
figure, where there are three segments for inclination i — 
0°, 30° and 60°; each segment consists of two values of 
the turbulent velocity, Sv — 0.5 and 2.0 km s~^. Similar 
sets for — 10^^ (middle part) and Xf — 10^^° (left 
lower part) are clearly in a different part of the diagram 
and do not overlap with each other. The variations due to 
inclination and turbulent velocity are smaller than those 
associated with freeze-out. 

4. CONCLUSIONS 

We have shown that it should be possible to observe 
some far-infrared rotational lines of water with the HIFI 
instrument on Herschel that are produced in regions be- 
tween radii i? ~ 10 — 100 AU. The current scheduled 
observations of the lio — loi and In — Oqo transitions 
would provide information on the spatial distribution of 
water and on turbulence if the inclination of the disk 
is known. The ratios of these lines are quite similar 
for the different models (~ 0.9 — 1.0), which is due to 
the high opacities in both lines. The line shapes of the 
lio ~ loi transition for different water vapor residuals in 
the freeze-out zone are quite similar at the same inclina- 
tion angle. When measurements of the lio — loi line are 
combined with higher transitions, such as the 3i2 — 3o3, 
and 221 — 2i2 lines, the abundance of the residual water 
vapor in the freeze-out zone can be determined (see Fig. 
2]). These Herschel observations of rotational lines will 
complement observations of shorter-wavelength excita- 
tion hues by Spitzer (10-20^m, e.g., Carr & Najita ,20081 : 
Salyk et al. l2008f ). SOFI A (6Mm), and near-infrared 
ro- vib rational transitions (jCarr et al.l 120041 : ISalvk et al.l 
120081 )). These lines trace the inner regions of the disk 
out to radii i? ^ 1 — 2 AU, with the potential of pro- 
viding information on the conditions in the immediate 
neighborhood of the snow line. 
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